Isofrequency pairing of geodesic orbits in Kerr geometry 
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Bound geodesic orbits around a Kerr black hole can be parametrized by three constants of the motion: the 
(specific) orbital energy, angular momentum and Carter constant. Generically, each orbit also has associated 
with it three frequencies, related to the radial, longitudinal and (mean) azimuthal motions. Here we note the 
curious fact that these two ways of characterizing bound geodesies are not in a one-to-one correspondence. 
While the former uniquely specifies an orbit up to initial conditions, the latter does not: there is a (strong-field) 
region of the parameter space in which pairs of physically distinct orbits can have the same three frequencies. 
In each such isofrequency pair the two orbits exhibit the same rate of periastron precession and the same rate of 
Lense-Thirring precession of the orbital plane, and (in a certain sense) they remain "synchronized" in phase. 



I. INTRODUCTION 

The motion of test bodies in the Kerr metric of a rotating 
black hole has been studied for almost half a century (see, 
e.g., 1 1-5 1). Much of the more recent work is motivated by the 
need to understand radiative inspirals into a Kerr black hole as 
sources of gravitational waves for future detector experiments. 
Examples of recent work include an action-angle formalism 
[|6 1, a frequency-domain method for computing functionals of 
the orbit (such as the gravitational perturbation from an orbit- 
ing test particle) Q, a system for classifying Kerr orbits |l8]|9l, 
and an analytic method for solving the geodesic equations of 
motion [lOJ. 

Timelike geodesies of the Kerr geometry are completely 
integrable. They admit three non trivial constants of motion 
("first integrals"), each associated with a Killing field of the 
Kerr background: the time-translation and rotational Killing 
vectors give rise to conserved (specific) orbital energy £ and 
azimuthal angular momentum £, and the second-rank Killing 
tensor discovered by Carter IT] gives rise to what is known as 
the Carter constant, Q. Up to initial conditions, these three 
constants of motion uniquely label all timelike geodesies of 
the Kerr geometry. 

This paper is concerned with the family of bound geodesic 
orbits. Each bound orbit is confined to the interior of a com- 
pact spatial torus given by rp < r < ra and 6'niin < ^' < 
TT — ^min, whcrc hereafter t, r, 9, (p are Boyer-Lindquist (BL) 
coordinates, r = rp, ra are two radial turning points ("perias- 
tron" and "apastron", respectively), and 9 — 0min,7r ^ 6*111111 
are two longitudinal turning points. Generically, the motion 
is ergodic ("space filling"), in the sense that any point inside 
the torus is visited by the orbit within a finite time t (excep- 
tional are "resonant" orbits, mentioned briefly below). The 
triplet {rp, ra, ^niin} provides an alternative parametrization 
of bound geodesies, which is in a one-to-one correspondence 
with that of {£■,£, Q}. 

Generically, bound orbits are triperiodic, with three fre- 
quencies n,r, and £7^ associated with the motions in the 
radial, longitudinal and azimuthal directions, respectively. Of 
these, ilr and ftg are "libration"-type frequencies, defined 
from the (average) radial and longitudinal periods, while ft^ 



is a "rotation"-type frequency, describing the average rate at 
which the BL azimuthal phase (p accumulates in time. We de- 
fine the above frequencies with respect to BL time t; this is 
useful for many purposes, because t is also the proper time 
of an asymptotically far static observer (e.g., a gravitational- 
wave detector). It is important to note that, in general, the 
orbital radius r and polar angle 9 of a given orbit are not (sep- 
arately) periodic functions of t: the t-interval between suc- 
cessive periastron passages is not constant, and the t-interval 
between successive 9 — ^i^in passages is not constant ei- 
ther. There is a choice of a time variable (the so-called "Mino 
time" — see below) in terms of which the radial and longitu- 
dinal motions completely separate and become precisely peri- 
odic. However, in terms of BL time t, the orbital periodicity 
can generally only be defined through an infinite time average 
(or, equivalently, through an average over the orbital torusQ). 
We shall define the BL-time frequencies more precisely be- 
low, following Schmidt |6| and Drasco and Hughes |7| 

The above general description simplifies in several special 
cases. If the ratio fir /fie is a rational number ("resonant or- 
bits"), then the trajectory traced by the orbit in the r-9 plane 
is closed (with a finite t-period), and the ergodicity property 
is lost. If the orbit is equatorial (6 — const = 7r/2), then fig 
loses its meaning, and the orbit becomes biperiodic with fre- 
quencies fir and n^p; in this case the radial motion is strictly 
periodic, with a radial period 27r/ilr. Similarly, if the orbit 
is circular (rp — ra), then fir loses its meaning, the orbit be- 
comes biperiodic with frequencies fig and fl^, and the longi- 
tudinal motion is strictly periodic with period 2Tt /fig. (Orbits 
that are both equatorial and circular are singly periodic with 
frequency fl^.) Finally, in the special case of a Schwarzschild 
black hole, one can always set up the BL system so that the 
orbit is equatorial and biperiodic with frequencies fir and fl^. 

The purpose of this article is to challenge the commonly 
held notion (see, e.g., ifTTI ) that the trio of frequencies 
{fl^, fir, fig} provides a good parametrization of generic 
bound geodesies in Kerr, one which is in a one-to-one cor- 
respondence with {£,£,2} or {rp, ra, ^min}- We show that 
this is not the case: there are infinitely many pairs of "isofre- 
quency" orbits, which are physically distinct (i.e., have dif- 
ferent {£,£,2} values) and yet they share the same values 
of {flip, fir, fig}- This point was already made briefly by two 
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of us in Appendix A of flT] in reference to a Schwarzschild 
black hole (where orbits are biperiodic, and two isofrequency 
orbits share the same values of il^p and O^). Here we first 
revisit the Schwarzschild problem to provide a further illumi- 
nation of this phenomenon, and then extend the analysis to 
the Kerr case, showing that isofrequency pairing occurs even 
among triperiodic orbits. 

We shall on occasion refer to a pair of isofrequency orbits as 
"synchronous", because the phases of such orbits remain syn- 
chronized in an average sense. For example, two equatorial 
isofrequency orbits that pass through their periastra simulta- 
neously at (f — will reach their next periastra at the same 
time and with the same azimuthal phase; they will have expe- 
rienced an identical amount of periastron advance. Although 
such orbits go "in and out of phase" between periastron pas- 
sages, their phase remains synchronized "on average". We 
will present some graphics to illustrate this behavior. 

Throughout this article we use geometric units such that the 
gravitational constant and the speed of light are both equal to 
unity. We denote the black hole's mass and spin by M and 
aM, respectively. We use an over-tilde to denote adimension- 
alization using ill; for example, (l^ :— Mfii^ and a :— a/M. 
We adopt a convention whereby a > and a < corre- 
spond to prograde and retrograde orbits, respectively, with C 
always positive. We use the term "orbit" synonymously with 
"timelike geodesic orbit". In Sec.[Il]we consider (biperiodic) 
synchronous orbits in Schwarzschild geometry (a = 0). We 
delineate the region in the parameter space where such or- 
bits occur, and also provide an intuitive explanation as to why 



one-to-one correspondence with that of {£, C}. Explicitly, 



isofrequency pairing must occur. In Sec. Ill we generalize our 
discussion to the Kerr case, where we consider first equatorial 
orbits and then generic, triperiodic orbits. 
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which can be inverted (for real e, p) to give unique expressions 
foip{£, C) and e(£, C). In the (p, e) space, bound orbits span 
the range < e < 1 withp > Ps{e) := 6 + 2e. The boundary 
Ps (e) ("separatrix") separates between stable and unstable or- 
bits in the (p, e) space |14|. The (p, e) = (6,0) terminus of 
the separatrix curve is known as the innermost stable circular 
orbit (ISCO). The existence of a separatrix is one of the salient 
features of motion in black hole spacetimes, and it marks a 
major qualitative departure from Newtonian dynamics. As we 
shall see, the occurrence of isofrequency pairing of orbits is 
intimately related to the existence of a separatrix. 

The function r{t) is periodic with (t-)period T^. Follow- 
ing Darwin 1 13 1, it is convenient to introduce the "relativistic 
anomaly" parameter %, which is related to t via 

dt _ Afp2[(p_2)2-4e2]i/2(p_6-2ecosx)"^/^ 
dx (p — 2 — 2ecosx)(l + ecosx)^ ' 

(4) 

and in terms of which the radial motion is given simply by 
r (x) = Mp/ (1 + e cos x) (taking x = at a periastron pas- 
sage). The radial period can then be computed via 



Jo dx 
with associated radial frequency 

2n 



Tr 



(5) 
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II. ISOFREQUENCY ORBITS IN SCHWARZSCHILD 
GEOMETRY 



The azimuthal frequency of the orbit is defined as the aver- 
age of dip/dt (with respect to t) over a complete radial period: 



A. Orbital frequencies and separatrix 

The radial motion of geodesic test particles in the equatorial 
plane of a Schwarzschild black hole satisfies 



£2 



(1) 



where a dot denotes differentiation with respect to proper- 
time, and V{r; C) is an effective potential for the radial mo- 
tion. Bound orbits exist for C > 2\/iM with ^ < £ < I. 
For each , £} in this range, r^(r) has three real roots, and 
motion is allowed between the second largest and largest of 
these, which we label r-p and 7-a, respectively. A convenient 
alternative parametrization of bound orbits is provided by the 
pair of values {p, e} defined through 



Mp 



2rpra 



(2) 



which are relativistic generalizations of semi-latus rectum and 
eccentricity, respectively |13|. This parametrization is in a 
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^dt = — 
dt Tr ' 
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where A(p is the azimuthal phase accumulated over time in- 
terval Tr- The latter can be computed via 



2tt 



/p 



^/p-6-2ecosx 



where e:=p~p,{e) andK{x) := J^^^ d0{l~xsm^ d)-^/^ 
is the complete elliptic integral of the first kind. 

At the separatrix limit, e ^ 0+, both Aip and Tr diverge at 
a similar rate [see Eqs. ^ and ( [T0| below], so that fi^ — ^ 
while fl^ attains a finite value [= {M/r'^)^/^, corresponding 
to the frequency of the unstable circular orbit of radius rp = 
p/{l + e)]. This gives rise to the well known "zoom-whirl" 
behavior ifTsll : orbits with e ^ 1 can "whirl" around the black 
hole many times near the periastron before "zooming" back 
out towards the apastron. 



dx 



(8) 
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B. Isofrequency orbits 

As pointed out in Ref. lfT2l . the Jacobian matrix of the trans- 
formation {p, e) — > {rtr, il.ip) turns out to be singular along a 
certain curve in the parameter space, well outside the separa- 
trix. This indicates that the transformation is not bijective. To 
see this most clearly it is instructive to move to a new orbital 
parametrization given by the pair (17^, e). This reparametriza- 
tion is admissible because (i) as argued above, the original 
parametrization {p, e) is a good one, and (ii) as can be eas- 
ily checked, il^p is a monotonically decreasing function of p 
for any fixed e. Our argument now follows from examining 
the structure of the Qr = const contour lines in the (17^,6) 
plane, as shown in Figure [T] The key feature here is that 
some ilr ~ const contours have vertical tangents (the locus of 
which is shown by the dashed black line in the figure). Each 
of these contour lines is intersected twice by vertical lines 
just right of the vertical tangent. But vertical lines are also 
= const contours, and so the two intersections mark a pair 
of isofrequency orbits. (Any two such isofrequency orbits are 
clearly physically distinct: they have different eccentricities.) 

In Fig. [2] we show, superposed, the orbital trajectories of 
a sample pair of isofrequency orbits of rather different ec- 
centricities. The radial and azimuthal motions of these two 
orbits are plotted in Fig.|3] Since the rate of relativistic perias- 
tron advance depends only on the frequency ratio fl^/ftr, two 
isofrequency orbits will exhibit the same rate of advance. This 
means that their phase remains "synchronized" on average, a 
behavior illustrated in the figures. 

Before giving a more detailed analysis, let us remark on 
the practicalities of producing the contour map of Fig.[T] The 
relation flr{ilip, e) is not known analytically, so we resort to 
a numerical calculation: First, for a given e, we numerically 
invert the relation e) [Eq. (j7]|] to find p{flip, e). Then 

we use Eq. (|6]l to obtain f7r(p(f7ip, e), e). Much of the inter- 
esting portion of the parameter space for our purpose lies very 
near the separatrix, where it becomes numerically challenging 
to evaluate the divergent quantities T^, and their ratio in 
Eq. (j7]i. In this problematic domain we instead use the near- 
separatrix analytic expansions Iil4l 
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+ 0(6 lege). (10) 



4(l-e2)3/2 

Here, the integral /(e) := /^^(l + e cos x)^^ D{cos x), with 
3 + 2e — cos^ X 



-D(cosx) 



^[2(1- cos x) 



1/2 



2 + e(l — cos x) 
3 + e- ^(7e-3)(l + cosx), (H) 



is easily evaluated numerically. 

It is in fact not hard to demonstrate the existence of isofre- 
quency orbits without resorting to a numerical calculation as 




FIG. 1. The (n^,e) parameter space for bound geodesic orbits in 
Schwarzschild geometry. Bound orbits are confined to the region 
right of the curve marked separatrix. Thin (blue) curves are contour 
lines of constant Q.r. The marginal contour line Q,r = is shown as 
a thick (red) line. Q,r takes its greatest value at the point marked c, 
representing a (slightly perturbed) circular orbit of radius = 8A/. 
The dotted (black) line shows the curve along which the Jacobian 
matrix of the transformation (p, e) ($1^,^^^) becomes singular. 
The singular curve intersects the e = axis at 6, corresponding to a 
circular orbit of radius n = (39 + 0l5)M/8 ~ 6.3802i\f. Any 
vertical (Q,^ — const) line left of b intersects some Q,,. — const 
contours twice. Each pair of intersections identifies a pair of isofre- 
quency orbits; a sample pair is marked in the plot. Each and every 
orbit between the separatrix and the singular curve has an isofre- 
quency dual between the singular curve and the dashed (green) curve 
marked COD (for circular-orbit duals). The COD is the locus of all 
orbits dual to circular orbits with radii 6AI = ri < r < rb- 



above. The argument follows from a few simple observations, 
which we now describe. First, it is easily established that, in 
the {e,fl^) plane, the separatrix e = es{^^) is a curve of a 
positive slope as shown in Fig. [T| (noting that in the figure we 
have chosen the horizontal axis so that fl^a increases to the 



left). To see this, use Eq. (j?]) with (j9]l and ( 10 1 to derive the 
relation ilip{e) along the separatrix, and invert to obtain 



1 - 2n 



1 

2/3 



(12) 



where, recall, Ct^p :— Mil^p. This gives dts/tKl^ > in the 
relevant range < 17^ < 1/8. Next, examine the curve ilr = 
in the (e, Vl^) plane: It runs up along the separatrix, then 
proceeds horizontally along the line e — 1 (which represents 
orbits with Ta — > oo and hence Tr — > oo), and finally descends 
along the line il^p = (which represents weak-field orbits 
with Tp — > oo, for which both frequencies vanish). Hence, 
the ilr = contour is represented by the thick red line in Fig. 
[T] circumscribing the parameter space of bound orbits on 3 
sides. From continuity, it is now clear that a contour line of 
sufficiently small flj. must "bend backward" inside the wedge 
formed by the separatrix and the e = 1 line, so that it becomes 
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FIG. 2. Orbital trajectories in the equatorial plane of a Schwarzschild 
black hole, for a sample pair of isofrequency orbits. The mo- 
tion is anticlockwise, and the black hole is drawn to scale. Orbit 
1 (red, round markers) has parameters (pi,ei) — (6.255,0.05), 
and orbit 2 (blue, square markers) has parameters (^2,62) — 
(6.718788076,0.3522488173). Both share the same orbital fre- 
quencies, (ClrM^) ~ (0.01257801,0.06426083). The orbital pe- 
riod of both orbits is Tr — 499.535318A/ and each accumulates 
Aif ~ 32.100669 radians during that period. Both orbits start at 
their periastron marker '0' along the radial line PI. Each successive 
marker shows the orbital phase after a time period of n x rr/8, where 
n is the marker number. At T,./2 (marker 4) both orbits are synchro- 
nized again at their apastra along the line Al. When each test body 
has completed one orbit (marker 8) they are again synchronized at 
their periastra along the line P2. Both orbits have precessed by the 
same amount over their common radial period. 



r{t)/M 




FIG. 3. Evolution of r{t) and if{t) — fl^pt for the isofrequency pair 
shown in Fig. [2] Both radial and azimuthal motions are "phase- 
synchronized" on average. 



vertical at a point. The existence of isofrequency pairs follows 
immediately, as discussed above. 

Let us now delineate the region in the parameter space 
where isofrequency pairing occurs. In Fig. [T]we have indi- 
cated in a dotted black line the curve along which the transfor- 
mation {p, e) o {fir, flip) becomes singular. Each and every 
orbit left of this singular curve has an isofrequency dual right 
of the curve. In particular, each and every circular (e — 0) 
orbit on the open segment (i, b) has an isofrequency dual on 
the dashed green line marked as circular-orbit duals (COD). 
(Here we define the radial frequency of a circular orbit to be 
that of a slightly eccentric orbit, at the limit e — > 0.) Hence, 
each and every orbit between the separatrix and the singular 
curve has an isofrequency dual between the singular curve and 
the COD, and vice versa. We conclude that (i) all isofrequency 
pairs are confined to the region left of the COD, and (ii) every 
orbit left of the COD has an isofrequency dual. 

How "strong field" is the region left of the COD, where 
isofrequency pairing occurs? The isofrequency pair of lowest 
azimuthal frequency sits where the singular curve intersects 
the e = axis, at point b (refer again to Fig.[T]l. To calculate 
the value of fl^ at b, we analytically Taylor-expand the Jaco- 
bian determinant J \d{flr, flip)/d{p, e)| in e about e = 
(for fixed p). We find, to leading order. 



J(e ^ 0) 



9(4p2 - 39p + 86) 

"47Vf2p9/2(p_2)(p_6)3/2' 

of which the relevant root is 



P 



6.3802. 



(13) 



(14) 



This corresponds to a circular orbit of radius rf, ~ 6.3802Af 
and frequency fl^ (M/rg)-^ ~ 0.06205/M. Recall this 
is the lowest frequency of any isofrequency pair. The isofre- 
quency pair of highest frequency sits at the upper-left corner 
of the diagram in Fig. [ij it has Cl^ = 1/8. Hence, for a 
Schwarzschild black hole, the range of isofrequency pairing 
is given by 



0.06205 < < 0.125. 



(15) 



(For comparison, the ISCO frequency is fl^ — 6^"^/^ ~ 
0.068.) Evidently, the phenomenon is confined to the very 
strong-field regime of the Schwarzschild black hole. 

Finally, we note that all orbits in isofrequency pairs are 
strongly zoom-whirling. For example, the lowest-frequency 
isofrequency pair mentioned above (slightly perturbed circu- 
lar orbits of radii r — ^ r^) have Aip ~ 4.1 x 27r, i.e., they 
each complete more than 4 full revolutions in cp over a single 
radial period. This behavior is also manifest in the example 
shown in Fig.|2] 



III. ISOFREQUENCY ORBITS IN KERR GEOMETRY 

A. Equatorial orbits 

We consider first the case of equatorial orbits, in which 
the treatment is entirely analogous to that of orbits in 
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Schwarzschild spacetime. Equatorial orbits have Q=0, and 
are therefore parametrized by the pair {£,C} alone. As in 
the Schwarzschild case, bound equatorial orbits may instead 
be parametrized by the (BL coordinate values of the) turning 
points {ra, r-p}, or by a pair {p, e} defined from them as in 
Eq. (|2|. One can then write integral expressions analogous to 
Eqs.^ [with (jSj and Q] and ^ [with ^] for the radial and 
azimuthal frequencies of the motion; the dependence upon 
the black hole's spin a only enters via the expUcit form of 
the functions e, a) and c?i^/(ix(x;?'i e, a), which 

are significantly more complicated than their Schwarzschild 
(a = 0) reductions. The integral formulas for fi^ and fi^, 
for arbitrary spin, can be found in Sec. II. A of Ref. ifTSl . and 
an analytic formula for the separatrix curve, Ps{e), again for 
arbitrary spin, is given in Ref. |16|. We will not reproduce 
these expressions here given their complexity, and since we 
will be giving explicit formulas for generic orbits in the next 
subsection. 

One finds that our intuitive argument for the exis- 
tence of isofrequency orbits carries over directly from the 
Schwarzschild case to equatorial orbits in Kerr. Along the 
separatrix of the Kerr black hole, the function es{ilip) is 
most neatly expressed in terms of the periastron radius fp = 
{fl^^ — a)^/'^ (which, on the separatrix, corresponds to the 
radius of an unstable circular orbit of frequency n^,) lfT6ll : 



-fp + 6fp — 8af 



3a 



;;2 



2r„ 



(16) 




flr=0 

sing, curve 
COD 



Ma, 



FIG. 4. The (Jl^ , e) parameter space for bound equatorial geodesic 
orbits in Kerr geometry with a — 0.5AI. Compare with Fig. ^ 
The relevant features are as in the Schwarzschild case, and the ex- 
istence of isofrequency pairing below the COD is similarly evi- 
dent. We indicate a sample pair with {p, e) — (4.915656, 0.45) and 
(4.62288270,0.26313140), both having frequencies (0^,0^) = 
(0.112675037, 0.01291945). Labelled points on the horizontal axis 
correspond to circular orbits of radii (left to right) — 3.8994M 
(whirl radius of marginally bound marginally stable orbit; orbit 
of highest azimuthal frequency), ~ 4.2330i\f (ISCO), n ~ 
4.5039A/ (outermost orbit in an isofrequency pair), and Vc — 
5.7628A/ (orbit of highest radial frequency, Milr ~ 0.03312). 



^ Qmax 

is the whirl 



It can be easily checked that dcg/df^ < and df^/dO.^ < 
for all a and all il^p in the relevant range < 
leading, again, to deg/dVi^ > 0. [Here SI 
frequency of the marginally bound and marginally stable or- 
bit with £ = 1 (and e = 1), an expression for which will 
be given in Eq. ( [TS] ) below.] The pattern of the Q,r = const 
contour lines in the (e, Vl^) plane should therefore be quali- 
tatively as in Fig. [T] including the crucial feature that contour 
lines "curve back" inside the wedge formed by the separa- 
trix and the e — 1 line. It follows that isofrequency pairing 
should be a feature of equatorial orbits for any black hole spin 
a (and, in particular, we expect to see it in both prograde and 
retrograde orbits). 

Figure |4] shows an actual contour-line map, similar to that 
in Fig. [T] for the sample case a = 0.5M. The S7j. = const 
contours were computed numerically as in the Schwarzschild 
case, this time using the integral expressions from Ref. ifTSl . 
Near the separatrix we have used the asymptotic expressions 
also given in ifTSl . Evidently, the essential features are as in 
the Schwarzschild case. One again identifies a singular curve 
and a COD curve in the (e, Vl^) plane, so that for any orbit 
between the separatrix and the singular curve there exists a 
dual isofrequency orbit between the singular curve and the 
COD, and vice versa. The situation is qualitatively the same 
for other values of the spin and for retrograde orbits. 

Let us identify the frequency range SI""" (a) < Vl^p < 
r2™^'^(a) where isofrequency pairing occurs. The a ^ ver- 
sion of Eq. ( 1 3 I is too complicated to be solved analytically for 
p = rb (the radius of the outermost circular orbit belonging to 



an isofrequency pair) as we have done in the Schwarzschild 
case, so we resort to numerical solutions. Table|l]lists rjj val- 
ues for a sample of black hole spins. Once a numerical value 



for is at hand (for a given a), S^""" is obtained via 



1 



^ -3/2 , 



(17) 



where we have used the general relation between the fre- 
quency of a circular equatorial orbit and its BL radius [2]. 
The maximal value Sl™^^ corresponds to the whirl frequency 
of the marginally bound marginally stable orbit with e = 1 
(top left comer in Fig.[4]i. It is given by 



Qmax 



1 



(2 - a + 2^/^^)3/2 + , 



(18) 



The range rj^'"(a) < < n"'"''' 



(a) is illustared in Fig.|5] 



B. Triperiodic orbits: frequencies and separatrix 

We now turn to consider generic bound motion in Kerr ge- 
ometry. Nonequatorial orbits possess a third frequency, ilg, 
associated with the longitudinal motion. It is not immedi- 
ately obvious how the three fundamental frequencies can be 
computed in practice, since the radial and longitudinal mo- 
tions are coupled in the usual BL-coordinate representation 
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TABLE I. Numerical values for 7-5, the BL radius of the outer- 
most circular orbit belonging to an isofrequency pair (cf. Fig. |4j. 

1] marks the 
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The frequency fi™'" of this orbit [given in Eq. 
lower end of the frequency range where synchronous pairing oc- 
curs. For comparison, the second column displays the ISCO ra- 
dius risco (elsewhere in this paper denoted r^); it is given by |2] 
hsco = 3-^-Z2-sign(a)[(3-Zl)(3-^Zl-h2Z2)]'/^whereZl := 

:;2^i/3[(i+a)i/3 + (i_a)i/3)]a„^j^2 - (3a2 + z2)i/2 



1+(1- 



Numerical values are truncated at the 5th decimal place, rounding up. 



[see Eqs. ([19]) and pO| below]. Schmidt ||6l was able to de- 
rive formal expressions for the fundamental frequencies us- 
ing angle-action variables in the Hamilton-Jacobi formalism, 
which circumvented the problem of coupling. Mino ifTTll ob- 
served that the radial and longitudinal motions can in fact be 
decoupled using a simple transformation of the time coordi- 
nate, and Fujita and Hikida HOI (building on work by Drasco 
and Hughes |f7l|) used this to obtain closed-form analytic for- 
mulas for the three frequencies. We give their formulas below 
in a slightly modified form. (Fujita and Hikida considered the 
cases \a\ 7^ M and \a\ — M separately. For brevity we repro- 
duce here only the nonextremal case; expressions for \a\ — M 
can be found in Appendix B of 1 10|.) 

To establish some necessary notation, let us begin with the 
r and 6 components of the geodesic equation of motion. For 
bound (£ < 1), nonequatorial ((?„iin 7^ 7r/2) orbits around a 
rotating (a 7^ 0) black hole, these can be written in the form 



2^2 



7(ri - r)(r - r2)(r - r3)(r - ri), 



■5^2 -2 2 / 2 2\l 2 



(19) 



(20) 



where z :=cos6', S :=r^+a^z^, 7 := 1— and an overdot 
denotes differentiation with respect to proper time along the 
geodesic. The roots of the quartic expressions on the right- 
hand sides are certain functions of Q; the radial roots 
are ordered as ri > r2 > > r4, and the roots ±2:_, ±2:4. 
satisfy |z_| < 1 and |z+| > 1. Bound orbits have rp = r2 < 
r < r\ = and \z\ < z_ = cos^min (the latter inequality 
corresponds to 6',ni,j < < tt ~ 0„iin)- We may introduce 
the parametrization {p, e, (?min}, where p, e are defined from 



FIG. 5. Range of isofrequency pairing (shaded area), as a function 
of the black hole spin, for orbits in the equatorial plane. The large- 
spin portion of the plot is shown separately in an inset for clarity. 
Isofrequency orbits are confined to the strong-field frequency regime 
> The frequency 17™^^ is the highest attainable by any 

bound orbit (at given M, a), corresponding to the whirl frequency of 
the marginally bound marginally stable orbit with £ = 1. 



Tp, Ta as in Eq. Q. The above roots are then most succinctly 
expressed (using a "mixed" parametrization) as 



COS! 



ri = Ta 



2+ = 1 



Mp 



a^-f SIT? 6',ni„ 



1/2 



l-e 



r2 = r-p 



Mp 



+ y/a^ - 4/3 



I 



, (21) 



(22) 



(23) 



where a := 2M — [r^ + r^) and (3 :— a^Q/ (jri^rp). 
Note that the r and 6 motions are coupled, due to the factor 



(r, 9) on the left-hand sides of Eqs. ( 19 1 and ( 20 1. This can 



be easily rectified by introducing a new time parameter A (of- 
ten referred to as "Mino time" in recent literature), satisfying 
A = S^^. In terms of A, the r and 6 motions decouple, and 
each becomes manifestly periodic, with A-frequencies and 
Tg, respectively. One can also define the azimuthal frequency 
as the average of dip/dX with respect to A, where in gen- 
eral the average needs to be taken over an infinite time. The 
three A-frequencies are given explicitly (for |a| 7^ M) by ifTOl 



2K{kr) 

1 /2 

7ra7 ' z+ 

£n(z_,fc£/) a 



X 1 



K{kg) r+ - r 



2M£r+ - aC 



rs - r+ 



(24) 
(25) 
(26) 



7 



where U{x,y) := J^^"^ d9{l - a:: sin^ - ysin^ 

is the complete elHptic integral of the third kind, r± := M ± 
i\f 2 _ (i2^ (jjg arguments of the elliptic functions are 



1/2 



''3 rp 



(27) 



and hereafter we use (+ ^ ) to denote a term formed by 
interchanging the + and — subscripts in the previous terms 
within the enclosing brackets. In Eq. (|26]l we have also intro- 



duced 

for A = {r, +, — }, with 

^ ^ (?-a-rp)(r3 -r±) 

(^■a - 7-3) (^-p - f±) 



n(/t^, fcr) 

K{kr) 



hr 



(28) 



(29) 



Finally, the t-frequencies are obtained from the A- 
frequencies via Q 



r ' 



T 



(30) 



where F is the average of dt/dX with respect to A. The latter 
is given explicitly (for |a| 7^ M) by ifTOl 



F = Aivrs 



£Q{1-Ge) , £ 



9 

7zl 



+ 77 ['"3(^'a + rp + rs) - r^rp + [r^ + rp + ra + ri)Fr + [r^ - r^){rp ~ ri)Qr 



2M 



{UP£ - aC)r+ - 2Ma^£ 

»'3 - r+ 



(31) 



where we have also introduced 



Qb := 



E{kB) 

Kike) 



(32) 



for B = {r, 0}, with E{x) J^^'^ de{l - x sin^ 6)^/^ being 
the complete elliptic integral of the second kind. 

Equation (jSO] ), with the necessary substitutions from Eqs. 
(21 i-(p9l), and ( [32) i, gives closed-form expressions for 
the fundamental frequencies il.r, f^e and fl^, given the pa- 
rameters {£, £, Q} as well as the corresponding parameters 
{p, e, 6'niin}- To complete the formulation, one requires a link 
between the two sets of parameters. Explicit expressions for 
{£,£, Q} in terms of {p, e, (?min} were derived by Schmidt 
in Appendix B of Ref. f6 1 (they are reproduced in a somewhat 
more concise form in Appendix A of ifTSll ). With this link, Eq. 



(30 1 can be used to compute the fundamental frequencies for 
a geodesic with given {p, e, ^min}- 

The separatrix between stable and unstable orbits is given 
by the condition rp — r^, which identifies the point where the 



inner turning point of the bound orbit is lost [recall Eq. (19 1]. 
It can be checked that this condition coincides with ilj. = 0, 
as expected (note k^. = 1 — and Fa — — Qr along the 
separatrix). Using Eqs. (|22li and (|23ll, with the link between 



{£, £, Q} and {p, e, (?min} from II6U18I. the condition rp = 
ra translated to a relation between p, e and ^niin, which can 
be solved numerically for p to obtain the separatrix surface 
p = Ps{e, 6'min)- We checked, using numerical examples, that 
this procedure for identifying the separatrix is consistent with 
the analytical method of Ref. |16| for equatorial orbits, and 
with the alternative numerical method of Sundararajan liT9l 
for generic orbits. 



C. Isofrequency pairing in triperiodic orbits 

We now seek to demonstrate the existence of isofrequency 
pairs of triperiodic orbits, i.e., ones sharing all three funda- 
mental frequencies {fir, ^e, ^^p}- Here our analysis will not 
be as complete as it was for biperiodic orbits. Rather, we will 
content ourselves with demonstrating by way of numerical ex- 
ample that such pairing does indeed occur. 

To this end it will suffice to inspect the contour map of 
51r=const curves in the (e, fig) plane, for some fixed value 
of ilip. For this, we need to be able to compute fir given 
{fie, e, fl^}. To achieve this in practice we take the follow- 
ing steps. First, we numerically invert, for given e, 6',nin, the 
equation fl^p (p) = const (in the example presented below we 
take the constant to be 0.14A/^^). For this we use a bisec- 
tion method, taking as initial guess the value p = Ps{e, ^min) 
obtained using the method described above. Once we have 
the trio {p, e, 6'min}, we calculate the corresponding values of 
ilg and rtj. using the analytic expressions presented above. 
We repeat these two steps for a great many values of e and 
(?min, making sure to achieve a good coverage of the param- 
eter space, particularly near the separatrix. The outcome of 
this procedure is a list of {fig, e, £7^} values for many orbits, 
all with our fixed value of il^. This dataset can then be used to 
create a contour map of f2r=const curves in the (e, fig) plane. 

We remark that the analytic formulation by Fujita and 
Hikida proves extremely useful for our purpose, because 
is can be readily implemented on a computer algebra plat- 
form such as Mathematica, which allows for high precision 
floating-point arithmetic. In our procedure, such high preci- 
sion is crucial near the separatrix, and it avoids the need to use 
asymptotic expansions as in Eqs. (j9]l and ( 10 1. 
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FIG. 6. Illustration of isofrequency pairing in triperiodic orbits. Thin 
solid (blue) lines are contours of constant Sir in the (e, 0,$) plane, 
for inclined eccentric orbits with fixed 17^ — 0.14. Here a = 0.7M. 
(The "empty" lower-right corner of the diagram lies outside the pa- 
rameter space of bound orbits.) Dashed vertical lines are sample 
r2e=const contours, along each of which we indicate a pair of isofre- 
quency orbits. The parameters of these three pairs are given in Table 
[ll| (sample pairs T, '2' and '3', from right to left). All essential fea- 
tures are as in Figs. JT] and |4] Similar contour maps can be obtained 
for other values of Qu, and a. 



An example with fl^p = 0.14 is shown in Fig. |6j We ob- 
serve that the essential features of the contour map are just as 
in Figs. [T] and |4] In particular, there are vertical (ile=const) 
lines that cross single J7r=const contours twice. Each pair of 
intersections represents a pair of isofrequency orbits sharing 
all three frequencies {fir, He, flip}. The existence of isofre- 
quency pairing in tiiperiodic orbits is thus established. Con- 
tinuity suggests that there should be a certain volume in the 
3-dimensional parameter space where isofrequency orbits re- 
side, but here we will not endeavour to identify the boundaries 
of this volume. 

We have indicated in Fig. |6] three sample pairs of triperi- 
odic isofrequency orbits, whose parameters we give in Table 
[ill The third sample pair is visualized in real space in Fig. |7] 
and in Fig.lslwe illustrate the synchronized evolutions of r{t), 
9{t) and Lp(t) for this pair. We note that in the case of triperi- 
odic orbits the "synchronization" is not exact (because the ?- 
and 6 motions are not separately periodic). Rather, the isofre- 
quency orbits are synchronized only in a long-time average 
sense. For example, if the two orbits pass their respective pe- 
riastra at i = 0, they may pass subsequent periastra at slightly 
different times, but the discrepancy should average to zero 
over a long time. One way to identify such behavior is by in- 
specting the difference between the orbital phases of the two 
orbits: the difference will remain qusi-periodic only if the two 
are isofrequency. The lower panel of Fig. [8] exemplifies this 
for the 9 phase. 

Before concluding, let us comment on the validity of our 
numerical algorithm, which, as already mentioned, involves 
delicate high precision computation of the orbital frequencies. 
To establish confidence in our results we tested our code in 
a number of ways. First, we checked that our code repro- 
duces all the double-precision-accurate results for {£, £, Q] 




¥IG. 7. A sample pair of triperiodic isofrequency orbits for a = 
0.7M. The orbits depicted correspond to 'sample pair 3' from Ta- 
ble |Il] (also leftmost pair in Fig. [6]l, with 'orbit 1' shown on the 
left and 'orbit 2' shown on the right. The top row shows the mo- 
tion in the (x, i/)-plane and the bottom row shows the motion in the 
(a;, 2:)-plane, where x = r cos sin 6/A/, y — r simp sin 9/ M and 
z — r cos6/M. The black hole is shown to scale. In both orbits the 
motion begins at f = = A at periastron, with ip = and 8 — -k 12. 
In integrating the geodesic equations we used the method of Drasco 
and Hughes [7], which avoids numerical difficulties near the orbital 
turning points. We show the portion of the orbits between A = and 
A = 30A/-\ 



(given {p, e, 6',nin}) tabulated in Ref. ifTSl . We also verified, to 
over one hundred significant figures, that the results of Fujita 
and Hikida's orbital frequency formulas (in the form given 
above) agree with the results of Schmidt's less explicit for- 
mulas |l6l. We further validated our equations using a direct 
numerical integration of the A-time geodesic equations in a 
few test cases. We were able to reproduce the analytically 
calculated A-frequencies and Tg to within 25 significant 
figures. (The quantities and F involve infinite time aver- 
ages and are therefore less easily tested in this manner.) 



IV. CONCLUDING REMARKS 

In this article we have shown that the three fundamental 
frequencies of bound geodesies in Kerr geometry do not con- 
stitute a good parametrization of the orbits in the strong-field 
regime. We identified a mapping between pairs of physically 
distinct orbits that possess the same set of orbital frequen- 
cies. A pair of isofrequency orbits are "synchronous" in that 
they exhibit the same periastron and Lense-ThiiTing preces- 
sion rates. All orbits in isofrequency pairs are confined to the 
very strong-field regime near the innermost stable orbit — cf. 
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orbit 1 orbit 2 


fundamental frequencies 


sample pair 1 


P 

e 

^min 


3.615857065600587178089905 3.4855158540000000000000 
0.068206017767752935160626 0.01201000000000000000000 
1.953894865146840010777339 1.8378366992075975844562 


Qr = 0.009040307723329 
Qe = 0.1162029753375 
= 0.1400000000000 


sample pair 2 


P 

e 

^min 


3.572207717388546694585166 4.3523435765502772064368261758 
0.0388932801825514054684027 0.27300090000000000000000000000 
1.943740959072074359824863 2.3444136677413832042020276247 


Qr = 0.006051252001160 
Qe = 0.1162584817374 
= 0.1400000000000 


sample pair 3 


P 

e 

^min 


3.80671950837597698109947211 4.477551959004760003175297459526 
0. 105336584613486946768869507 0.300000000000000000000000000000 
2.11092907046831122994268532 2.395463898362217344327765579750 


Qr = 0.005364669707792 
Qe = 0.1163492371285 
= 0.1400000000000 



TABLE II. Sample pairs of triperiodic isofrequency orbits for a — 0.7M (cf. Fig.|6j. The high precision of the parameter values presented is 
necessary for the orbital frequencies to match to within the 13 significant figures displayed. This level of precision is required because, near 
the separatrix, small changes in the orbital parameters can result in comparatively large changes in the frequencies. 
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t/M 



cos[0i(t)] 
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FIG. 8. Evolution of r(t), ip(t) ~ Q.^t and cos[6li(t)] - cos[6l2(t)] 
for 'sample pair 3' of Table [HI and Fig.|7] Both orbits begin at f = 
at periastron with ip = Q and 9 — n/2. Triperiodic isofrequency or- 
bits are "synchronized" only in a long-time average sense. Periastra 
are reached only approximately at the same time (as a closer inspec- 
tion of the upper panel would reveal) but the time differences should 
average to zero over a long time. The same applies to the average az- 
imuthal motion (middle panel, where a close inspection reveals that 
the azimuthal phases of the two orbits are not in precise agreement at 
the periastra), and to the motion in 6 (lower panel). In the latter case 
we show the difference between the two longitudinal phases, which 
remains quasi-periodic. It would have not remained quasi-periodic 
had the two orbits not been in an isofrequency pair. 



Table [l] and Fig.|5] (Some orbits in isofrequency pairs have 
very large eccentricities and apastra at arbitrarily large radii, 
but their periastra are in the very strong field.) Our numerical 
experiments suggest that all members of isofrequency pairs 



are of "zoom-whirl" type, but this is yet to be checked more 
thoroughly in the case of triperiodic orbits and across all spin 
values. 

The first practical lesson from our analysis is a caution- 
ary note for colleagues studying the data-analysis problem 
for gravitational-wave detectors, in particular the problem of 
parameter extraction for systems of extreme-mass-ratio inspi- 
rals (EMRIs). The fundamental frequencies extracted from a 
"snapshot" of an EMRI waveform, on their own, as a mat- 
ter of principle, do not necessarily provide enough informa- 
tion from which to extract the system's intrinsic physical pa- 
rameters £, £, Q (or p, e, 6'min)- If the system is sufficiently 
close to the innermost stable orbit, a measurement of the in- 
stantaneous frequencies could at most narrow down on two 
possible sets of system parameters. This "degeneracy", how- 
ever, can be removed in any one of the following ways: (i) by 
examining the power spectrum of the waveform (the power 
distribution among the various harmonics of the fundamental 
frequencies will be different for the two orbits); (ii) by in- 
specting the waveform snippet in the time domain (the shape 
of the waveform is strongly dependent upon the eccentric- 
ity, for instance); or (iii) by accounting for radiation-reaction 
evolution effects (two orbits which are instantaneously isofre- 
quency will evolve radiatively in different ways). 

At a more fundamental level, our analysis identifies a new 
feature in the strong-field dynamics of compact-object bina- 
ries in general relativity. The fundamental frequencies in 
a bound binary (of any mass ratio) are important invariant 
characteristics of the "conservative" sector of the dynamics. 
As such they have long been studied in the context of post- 
Newtonian (PN) theory. The instantaneous frequencies in a 
binary of inspiralling black holes can even, nowadays, be ex- 
tracted from high-precision fully nonlinear simulations in nu- 
merical relativity (NR) — see, for example, Ref. f20l. Our 
analysis here revealed the occurrence of isofrequency pairing 
in the test-particle limit (i.e., the limit of vanishing mass ratio), 
but it is not unreasonable to speculate that the phenomenon 
is a general feature of the dynamics in strongly gravitating 
binaries, and would reveal itself also when the mass ratio is 
finite. It is not clear if available PN theory can predict isofre- 
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quency pairing — this would be interesting to check. When 
new, higher-order PN terms are calculated in the future, it 
would again be interesting to check if they reveal the phe- 
nomenon, as a way of assessing the faithfulness of the PN 
expressions in the strong-field regime. It would also be in- 
teresting to examine whether the phenomenon manifests itself 
in NR simulations of inspiralling black holes of comparable 
masses near the innermost stable orbit. 

Because the fundamental frequencies are invariant charac- 
teristics of the conservative dynamics, they are useful as ref- 
erence points for comparing between the predictions of dif- 
ferent approaches to the relativistic two-body problem. Re- 
cent examples of such "cross-cultural" comparisons include 

(i) calculations of the ISCO frequency in the self-force (SF), 
PN and effective-one-body (EOB) approaches |2T-53l; and 

(ii) calculations of the periastron advance in slightly eccentric 
orbits in SF, PN, EOB and NR CO] |22l El. hi both exam- 
ples (which involve two nonrotating black holes) relations be- 
tween the two invariant frequencies associated with infinites- 
imally perturbed circular orbits were utilized as benchmarks 
for comparison. The singular curve/surface identified in our 
current work is an invariant structure in the parameter space, 
which provides yet another, independent, comparison point in 
the strong field, this time utilizing eccentric orbits. 

As a first example, one could consider the function 
along the singular curve in the parameter space of nonrotating 
binaries (Fig. [T] shows this curve in the test-particle limit). In 
principle, one could compute this function in the SF approx- 
imation (i.e., order by order in the mass ratio), and perhaps 
also in fully nonlinear NR, making for an interesting com- 
parison. There may be a way of using the results of such a 
calculation to calibrate the potentials of EOB theory in the 
strong field, although how this could be done in practice is yet 
unclear ll25l . Comparison with existing PN expressions could 
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test the performance of the PN expansion in the strong field. A 
more constructive synergy could be achieved within the recent 
"phenomenological" approach to PN calculations, whereby 
high-order terms in the PN expansion are determined by fit- 
ting to numerical data from SF or NR calculations [26. ,27]. A 
faithful phenomenological PN model would need to be able to 
recover the singular curve in the strong field, perhaps through 
the inclusion of suitable "poles" in PN expressions. 

Finally, let us mention the intriguing possibility that isofre- 
quency pairing in astrophysical black holes (e.g., between 
clumps of accreting matter) could have observational impli- 
cations. The question is worth asking because we are at an 
era where astronomical observations in a range of electro- 
magnetic wavelengths routinely peer into processes deep in 
the strong-field potentials of accreting black holes. Quasi- 
periodic oscillations (QPOs) in x-rays from accreting black- 
hole systems probe the innermost regions of accretion disks 
ll28l . and (to a lesser extent) so do x-ray flares from the galac- 
tic center |29 1. Could the peculiar strong-gravity phenomenon 
of isofrequency pairing have a dynamical effect on matter or- 
biting the black hole, perhaps through resonant interaction? 
Although admittedly far-fetched, this possibility deserves ex- 
ploration. 
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